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Relativistic emission lines from accreting black holes 
The effect of disk truncation on line profiles 
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Abstract. Relativistic emission lines generated by thin accretion disks around rotating black holes are an important 
diagnostic tool for testing gravity near the horizon. The iron K-line is of special importance for the interpretation 
of the X-ray emission of Seyfert galaxies, quasars and galactic X-ray binary systems. A generalized kinematic 
model is presented which includes radial drifts and non-Keplerian rotations for the line emitters. The resulting 
line profiles are obtained with an object-oriented ray tracer operating in the curved Kerr background metric. The 
general form of the Doppler factor is presented which includes all kinds of poloidal and toroidal motions near the 
horizon. The parameters of the model include the spin parameter, the inclination, the truncation and outer radius 
of the disk, velocity profiles for rotation and radial drift, the emissivity profile and a multi-species line-system. 
The red wing fiux is generally reduced when radial drift is included as compared to the pure Keplerian velocity 
field. All resulting emission line profiles can be classified as triangular, double-horned, double-peaked, bumpy 
and shoulder-like. Of particular interest are emission line profiles generated by truncated standard accretion disks 
(TSD). It is also shown that the emissivity law has a great influence on the profiles. The characteristic shoulder-like 
line profile observed for the Seyfert galaxy MCG-6-30-15 can be reproduced for suitable parameters. 
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1. Introduction 

Emission lines originating from regions where the influ- 
ence of curved space-time can not be neglected are found 
in accreting black hole systems like Active Galactic Nuclei 
( AGN) , microquasars and Galactic Black Hole Candidates 
(GBHC). These astrophysical systems exhibit hot plasma 
falling into the deep potential well of black holes. The 
ion species in the plasma still have some electrons on 
lower shells so that there are several transitions between 
electron states that generate emission lines in the X-ray 
range. The most prominent line is the Fe Ka line at 6.4 
keV rest frame energy (near- neutral iron). Additionally, 
the Fe K/3 line at 7.06 keV and rarely Ni Ka at 7.48 
keV and Cr Ka at 5.41 keV can be applied to fit X-ray 
emission line systems. The fluorescence process and its 
transition energies depend on the ionization state of the 
species. One can distinguish hot and cold regimes of line 
emi ssion which can be desc ribed by an ionization parame- 
ter ijBallantvne et al.l200l|) . Generally, X-ray astronomers 
observe a multi-component emission line complex of these 
species. The Fe Ka line dominates the line system due to 
its high relative strength and is often indicated in unre- 
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duced global X-ray spectra. 

Broad X-ray emission lines originate from the hot inner 
parts of accreting black hole systems, typically from radii 
comparable to the radius of marginal stability. Naturally, 
these lines hint for a proximity of hot and cold material: 
hot material can be found in an optically thin corona, a 
place where cold seed photons are Comptonized to hard 
radiation. It is expected that the corona becomes opti- 
cally thick at super-Eddington accretion rates and is dom- 
inated by Comptonization. This latter scenario seems to 
hold for narrow-line Seyfert-1 galaxies (NLSls), whereas 
classical broad-line Seyfert-1 galaxies (BLSls) accrete at 
sub-Eddington to Eddington rates. 

Relatively cold material hides in the radiatively-efhcient, 
geometrically thin and optically thick accretion disk, the 
so call ed standard disk, hereafte r SSD {Shakura-Sunyaev 
Disk) llShakura fc Sunvaevlll97.'^ . The resulting structure 
equations for this standar d disk have been generali zed to 
the relativistic case later ()Novikov fc Thornelll974|) . The 
hard coronal input radiation usually described by a power 
law illuminates the cold disk. As a consequence, the sur- 
face of the disk is ionized in a thin layer. In a fluores- 
cence mechanism the hard radiative input is absorbed 
(threshold at 7.1 keV) and re-emitted in softer fluores- 
cence photons. But in the dominant competitive process. 
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the auger effect, the excited ions emit auger electrons 
and enrich the environment with hot electrons. One of 
the first applications of this X-ray fl uorescence was per - 
formed at the microquasar Cyg X-1 l)Fabian et alJll989|) . 
These observations were supplemented by radio-quiet 
AGN, e.g. Seyfer t galaxies, showing the same X-ray refiec- 
tion signatures llPqunds et al.l Il990l: iTanaka et alJ llQQSt 
lR,evnolds fc Nowakll2n0.'^ " 

But there are also sources exhibiting narrow X-ray emis- 
sion lines. The narrow lines do not show distortion by 
relativistic effects. It has been proposed that narrow X- 
ray emission lines can be explained by r eflection at the 
large-scale dusty torus ijTurner et al.ll99^ . Here, the nar- 
row Hne is overimposed on broad line profiles. This inter- 
pretation has been applied t o some sources, lik e Seyfert 
galaxies, e.g. MCG-6-30-15 llWilms et al.ll2nniK MCG- 
5-23 -16 llDewangan et al.1 12003|) and Quasars, e.g. Mrk 
205 (|Reeves et al.ll200l|) . 

The absence of the line feature in several Sefert-1 galaxies 
( Pfefferkorn et al.lf2fl0ll) can be explained by global accre- 
tion models. The first possible interpretation is that the 
inner edge of the disk seems to be far away from the black 
hole. If the fluorescence process is still possible in that 
distance, the skewness of the X-ray emission line vanishes 
and becomes rather Newtonian-like. In even larger dis- 
tances, it is supposed that the emission Hne feature disap- 
pears completely because irradiation and fluorescence is 
reduced. Accretion theory states that this scenario holds 
for low accretion rate: then the transition radius between 
standard accretio n disk and optically thin disk is large 
The mass determinations for Seyfert 
black holes confirm this interpretation while being in the 
intermediate range of 10*^ to 10^ Mq. These lighter super- 
massive black holes seem to accrete very weakly. 
X-ray spectra of AGN and X-ray binaries are dominated 
by broad Comptonized continua. The properties of these 
spectra are discussed in the next section. The first step 
for astronomers is to reduce the observed data files and 
to extract the line profile by subtracting this global con- 
tinuum emission. This procedure is already complicated 
because the intrinsic curvature of the Comptonized radi- 
ation is still topic of an ongoing debate. 

The resulting emission line profile serves for diagnos- 
tics and to constrain essential parameters of the accreting 
black hole system. So, astronomers found an X-ray im- 
print of black holes and can derive properties concerning 
the relative position from the inner accretion disk {stan- 
dard disk) to the observer, the inclination i or the rota- 
tional state of the black hole, the Kerr rotational pa- 
rameter a, the extension of the geometrically thin and 
optically thick standard disk, like inner radius rjn and 
outer radius Tout, the velocity field of the plasma, 
separated in toroidal w*, poloidal and radial com- 
ponents and the relative position of hot corona to cold 
accretion disk. 

Let us follow up the theoretical steps from black hole the- 
ory needed here until one gets simulated emission line pro- 
files: Sect. 01 introduces the metric of rotating black holes, 
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Fig. 1. Prototype X-ray spectrum of a type-1 AGN (il- 
lustration idea taken from Fabian, 1998). The spectrum is 
dominated by the global Comptonized continuum. At low 
energies, there is the soft excess, in principle a composi- 
tion of some black bodies belonging to the disk divided 
into rings. At energies around 1 keV, one can recognize 
a complex of absorption dips, originating from the warm 
absorbers. The reflection component consists of a broad 
bump peaking around 20 keV and the prominent fluores- 
cence emission line complex around 6.4 keV. 



the Kerr geometry, and provides the standard co-ordinate 
system in Boyer-Lindquist form. Then, the generalized 
GR Doppler factor is presented. The derivation of this 
quantity and its signification in GR ray tracing is elabo- 
rated, too. 

In Sect. 01 realistic plasma kinematics resulting from hy- 
drodynamic and magnetohydrodynamic accretion theory 
is discussed. The Kerr geometry has special features that 
are valid independently from the accretion model! This is 
accessable by analytical investigations. 
Another basic ingredient to calculate emission lines and 
spectra in general is the emissivity law that determines 
the emission behaviour of the line emitting region. In Sect. 
O alternative emissivity shapes besides th classical single 
power law are applied. 

Ray racing provides in a first step disk images. It may 
be possible to resolve these images observationally in near 
future. In any case, disk images as presented in Sect. |S1 
are very instructive unrevealing GR effects and influence 
of plasma dynamics. 

In Sect.d the calculation and basic properties of the emis- 
sion lines are discussed. 

In the final Sect. (HI the parameter space of relativistic 
emission lines is investigated and suitable criteria to fix 
line features are presented. One can apply these criteria 
to observed line profiles and derive a classification in the 
line jungle. 
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Fig. 2. Illustration of topological elements in the inner- 
most region of accreting black hole systems. The size and 
morphology of these elements depend on accretion rate, 
black hole mass and radiative cooling. 

2. X-ray spectra and topology of type-1 AGN 

The data quality increased conspicuously with the launch 
of the space-based X-ray observatories Chandra and 
XMM-Newton. The challenge in interpreting X-ray spec- 
tra of AGN and XRBs is to correctly reduce the data and 
subtract the continuum. Afterwards one can discuss the 
formation of other spectral components, such as emission 
lines. 

X-ray spectra of Seyfert galaxies have been investi- 
gated for several years. These types of AGN were the 
first objects that ga ve evidence for rotating black holes 
i Tanaka et al.lll995[) and confirmed the paradigm of the 



AGN engine: an accreting supermassive black hole pro- 
ducing high multi- wavelength luminosity. Only a few years 
later, high-resolution observatories allowed the extraction 
of iron K emission li nes in more distant objects, Quasars 
l Yaaoob et alJlTflflfll^ . In the course of the unification 
scheme of AGN, the different AGN types are consid- 
ered as different evolution phases of a prototype AGN 
with some pecularities that mainly arise by differencies in 
orientation of the accretion disk, accretion rate and black 
hole mass. Especially, the discrimination between Seyferts 
and Quasars has become somewhat arbitrary. Only the 
amount of radiation, w 10^^ Lq for Seyferts and lO^'^ to 
10^^ Lq for Quasars, is the main difference between these 
two subclasses. Concerning the iron emission line feature, 
it can be stated that Quasars exhibit increased ionization 
of the standard disk surface due to enhanced central lu- 
minosity. Therefore, weaker emissio n line signatures are 
expected llR.evnolds NowaklEini. 
The dichotomy phenomenon caused by disk orientation 
and dust torus obscuring, a generalization that can be 
applied to Seyferts and Quasars, is well-known: one can 
distinguish type— 1 (low inclination) and type— 2 AGN 
(high inclination). Hence, it is for observational reasons 
only possible to study relativistic emission lines of AGN 
of type 1. 

X-ray spectra of Seyfert galaxies in the range from 0.1 to 
several 100 keV show many features as illustrated in Fig. 
n] overlaying the complete X-ray range there is a con- 



tinuum originating from Comptonized radiation. This 
is the direct radiation that reaches the observer coming 
from a hot corona {T « 30 keV). The photon number 
fiux per unit energy is usually modelled with a power 
law distribution, Ff^{E) cx E^^ . The soft input radia- 
tion from the cold accretion disk is reprocessed via un- 
saturated inverse Compton scattering by ultra-relativistic 
electrons in the corona. Fig. [21 depicts in a simple way 
one probable geometry of an accreting black hole systems 
that hold especially for AGN and microquasars, known 
as " sphere-l-disk geometry" . Matter has angular momen- 
tum and therefore spirales down in a fiat and cold stan- 
dard disk. The argument for the fiatness is because of 
vertical collapse due to efficient radiation cooling. The 
relati ye scale heigth H/r of the se standard disks is only 
10"^ ijShakura fc SunvaevlllQTl . The matter within this 
cold disk emits a multi-color blackbody spectrum. This 
is simply the superposition of a sequence of Planck spec- 
tra, e.g. rings each with a temperature T;. The innermost 
ring is the hottest one. Tin > ?! > Tout- The accretion 
rate determines the transition radius and this one fixes 
the temperature of the inner edge of the disk. At high ac- 
cretion rates, the inner edge is significantly hot so that a 
soft excess can be observed around 1 k eV. At sma ll ac- 
cretion rates the soft excess lacks (Esin et allllioTl) . The 
cold and thin accretion disk delivers the soft seed photons 
that are Comptonized in the hot, optically thin corona. 
The topology of the corona is still one of the essential 
open questions in X-ray spectroscopy. Different assump- 
tions such as slab and pa,tchv corona models have been 
made l R,evnolds fc Nowakl200,'^ . The observed hard spec- 
tra suggest rather "sphere-|-disk geometries" because slab 
and patchy corona models are more efficiently radiative- 
cooled by seed photons from the cold disk. The reverber- 
ation mapping technique may enlighten the spatial po- 
sition of corona to accretion disk. Theoretically, it will 
be task of upcoming relativistic MHD simulations to un- 
reveal the corona-disk geometry depending on accretion 
rate and radiative cooling. Apparently, it turns out that 
the coronal emissivity is signific antly enhanced at the in- 
ner edge of the standard disk ((Merloni fc Fabianl bOO,*^ 
and that MHD dissipation plays a crucial role to invalidate 
the zero-torque boundary condition. Then, magnetically- 
induced torques at radii comparable to the radius of 
marginal stability are expected to increase dissipation. 
This mechanism may provide the higher disk emissivities 
e(r) (X r~ *-^ as recently obse rved in Seyfert galaxy MCG- 
6-30-15 ijWilms et al.l l2C 
Here it is assumed that the corona is almost identi- 
cal with an advection-dominated torus that forms by 
feeding from the standard disk. An alternative may be 
the quasi-spherical rather Bondi- a.ccreting Advection- 
Dominated Accretion Flow (ADAF) ((Naravan fc YilllQQ^ 
that forms under other circumstances in accretion the- 
ory, e.g. depending on the accretion rate. The transi- 
tion region where the standard disk adjoins the hot torus 
is unstable. Hot corona and cold disk may constitute 
a sandwich configuration where the inner edge of the 
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cold disk oscillates in radial direction l)GraciJ200j) . This 
ADAF-SSD transition region seems to depend on the efh- 
ciency of different cooling channels as can be found in hy- 
drody namic simulations ijManmoto fc Katoll2000t iGracial 
horn . More recent pseudo-Newtonian hydrodynamic sim- 
ulations incorporating radiation cooling via synchrotron 
radiation, bremsstrahlung and Comptonization with the 
consideration of conductio n in a two-component plasrn a 
suggest disk truncation i Huieirat fc Camenzindll200(ll) . 
This means that the accretion disk can cut off even at 
radii larger than the radius of marginal stability, rmg. 
Truncation of standard disks (TSD) in general soft- 
ens the argument for observed evidences of rapidly spin- 
ning black holes. This is because truncated disks do not 
extend to inwardly. However the decreasing r^^ for 
faster spinning black holes was always the argum ent for 
evidence of Kerr black holes ijiwasawa et al.lll996l) . 
Particularly, the cooling time-scale of synchrotron radi- 
ation is typically in the millisecond domain and cools 
very fast the hot accretion flow. The Truncated Disk - 
Adve ctive Torus (TDAT) scenario l|Huieirat fc C amenzind' 
12000^) produces a hot inner torus that is not stable either. 
It is destablized by accretion and matter free-falls into the 
black hole. But a small fraction of matter can escape the 
black hole and is transported along magnetic field lines 
in polar regions. This drives outflows that consist of disk 
winds on one hand and of ergospheric Poynting fluxes on 
the other hand. This outflow is supposed to feed the large- 
scale jets of AGN or the blobs of microquasars. In a unified 
view, the total accretor mass, M, and the total mass ac- 
cretion rate, M, are strongly favoured to determine the 
spectral state and the jet injection rate and the continu- 
ity of the outflow. Relativistic emission lines may serve as 
a diagnostic tool to flx these parameters. 
The continuum of the hot corona reveals an intrinsical 
curvature and cuts off at the high-energy branch at sev- 
eral hundreds of keV. The cut-off a t SknTg indicates the 
plasm a temperature in the corona l|Rvbicki fc LightmanI 
At the low energy branch, some sources show a 
soft excess with a bumpy proflle. This can be identified 
with a multi-color blackbody originating from the cold 
Compton-thick accretion disk. At 1.0 keV there is a com- 
plex structure of warm absorbers, typically a combina- 
tion of adjacent absorption edges of a variety of species 
(CV, OVII, OVIII, NelX etc.). These dips and their 
potentially rel ativistic broadening are still t opic of an on- 
going debate ijMason et al.ll200,'jlLe3f200ll) . The reflec- 
tion component consists of a system of emission lines, 
most prominently the iron Ka at 6.4 keV typically, and 
a broad bump peaking at 20 to 30 keV. The reflection 
bump is a consequence of hard radiation from the corona 
hitting th e cold accretion disk that is then reflected to the 
observer l(R,evnoldjll9"9l . The origin of the iron Ka line 
is well understood: at typical temperatures of 10^ to 10^ 
keV in the inner accretion disk, iron is ionized but not 
completely stripped; the K- and L-shell are still popu- 
lated. First, iron is excitied by photo-electric absorption. 
The threshold lies at ~ 7.1 keV so that only hard X-ray 
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Fig. 3. Radial behaviour of the Boyer-Lindquist func- 
tions at extreme Kerr (a — 1.0) in the equatorial plane. 
Relativistic units were used. 

photons from the corona can release this primary process. 
Then, there are two competing processes: the dominant 
process (66 % probability) is the Auger effect where the 
excitation energy is emitted with an Auger electron. This 
mechanism is non-radiative and enriches the disk plasma 
with electrons. The second and here essential process is 
fluorescence (34 % probability): one electron of the L- 
shell undergoes the transition to the K-shell accompanied 
by the emission of the iron Ka fluorescence photon with 
6.4 keV rest frame energy^. The emission line energy and 
the existence of the fluorescence line in general depen d on 
the ionization state of the material ijRoss et alJll999t) . 
The significance of other elements contributing to the re- 
flection component is rather low, because the abundance, 
relative strength and fluorescenc e yield is low compared 
to that of iron jR,evnoldslll99i . At most, the iron Kff 
at 7.06 keV or nickel at 7.48 keV ijWang et all Il999(l 
and chromium at 5.41 keV could contribute marginally. 
Besides, the accretion disk has only a thin ionization layer. 
The line producing layer has only the depth of 0.1 to 1% of 
the disk thickness ijMatt et al.lll997(l . Therefore the con- 
tributions of ionized iron from deeper regions - or even a 
radiation transfer problem from lower ionized slabs - can 
be neglected due to high optical depths. 

3. Ray tracing in the Kerr geometry 

The Kerr metric describes rotating uncharged black holes. 
This vacuum solution of the Einstein fleld equations l|KerJ 
Il963l) generalizes the static Schwarzschild solution. The 
line element fuUfiUs the standard form of stationary and 
axisymmetric space-times 



ds^ 



„2* 



dt^ 



e2*(d$ - ujdt)^ ~ e^f'-dr^ - e^^^dO^ (1) 



The Kerr geometry in Boyer-Lindquist form is determined 
by the line element 



-a^dt^ + u^id^ - ujdtf + (p^/A) dr^ + p^dO^, (2) 



^ We neglect here any distinction of the Kai line at 6.404 
keV and Ka2 line at 6.391 keV as can be found in atomic 
physics. 
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with the canonical Boyer-Lindquist functions 



A 

= r 



- 2Mr + a^, 



2 ' a^cos^e. 



{r^ + a^f -a^ A silled, 

— sinfl, 
P 

2aMr 



S2 



(3) 

(4) 
(5) 
(6) 

(7) 
(8) 



M is the mass and a the spin parameter of the black hole. 
The Boyer-Lindquist co-ordinates are pseudo-spherical. 
The expressions above are called frame-dragging fre- 
quency or potential for angular momentum to, cylindrical 
radius lj, lapse function or redshift factor a and some ge- 
ometrical functions A and p. 

The rotation of the black hole can be parametrized by its 
specific angular momentum a, the so called Kerr param- 
eter. Relativistic units G ~ M — c ~ I arc used unless 
otherwise stated. Then, the natural length-scale is the 
gravitational radius, Tg = GM/c^, and a varies between 
— 1 and 1. The interval a € [—1,0] involves retrograde ro- 
tation between black hole and disk, the interval a £ [0, 1] 
means prograde rotation. From black hole theory, it is gen- 
erally possible that a takes all values in this interval, but 
accretion theory suggests that there is mainly one type of 
black holes: rotating Kerr black holes near its maximum 
rotation (a « 1). This is because black holes of mass M 
are spun up to a ~ Am/M by accreting a mass Am. At 
least the supermassive black holes in AGN that accrete for 
long times are expected to rotate very rapidly, a > 0.9. 
It is possible to extract ang ular momentum froni black 
holes via Blandford- Z na i ek llBlandford fc Znaiekl Eml) 
and Penrose processes l)Penrose fc Flovdlll97lll . But these 
effects may not be efficient enough to slow down the black 
hole rotation significantly. 

Fig. El illustrates the radial dependence of the Boyer- 
Lindquist functions in the equatorial plane {6 = 7r/2) at 
the maximum value of black hole rotation, a = 1. One can 
see the frame-dragging frequency lu that increases signif- 
icantly at smaller radii and reaches oj{rii) = fin = 1/2 
in relativistic units at the horizon. This behaviour shows 
that dragging of inertial frames is important when reach- 
ing the ergosphere. Another essential feature is the steep 
decrease of the redshift factor a. This is mainly the cause 
that emission is strongly suppressed when the distance be- 
comes smaller than approximatly the radius of marginal 
stability, rms, which satisfies 



M [3 + Z2T y/{3-Zi){3 + Zi + 2Z2] 

1/3 



(9) 
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Fig. 4. Schematical representation of a Kerr ray tracer. 
The light rays start at the screen {back tracking) and hit 
the disk surface in the equatorial plane. On the screen the 
lensed image is formed as seen by a distant observer. Our 
solver folds directly from equatorial plane to screen. 



(upper sign prograde, lower sign retrograde orbits). 
Our obi ect-oriented c ode KBHRT {Kerr Black Hole Ray 
Tracer) ^lMiilleI^l20n(]^ programmed in C-I--I- is used to de- 
rive in the first step the image of a relativistic disk and 
in the second step to calculate the line flux by integrating 
over this image and weighting with the relativistic general- 
ized Doppler factor and the radial disk emissivity profile. 
The imaging of the disk is not observationally possible 
because one can not (yet) resolve these distant accreting 
black hole systems. But these images show nicely the in- 
fluence of curved space-time on photon paths. The flux is 
indeed the observationally accessible information and can 
be determined afterwards. 

The code which is based on a C-code l|Fanton et alJ 
I1997II calculates the photon path by integrating the GR 
geodesies equation using the complete set of integrals 
of motion. Besides the canonical conservative quanti- 
ties like mass, energy and angular momentum, the Kerr 
geometry exh ibits a fourth one: the Carter constant 
l|Carteil 1968|l . This quantity is associated with the ra- 
dial and poloidal photon momentum. The second order 
null geodesies equation reduces by means of these four 
conservatives to a set of four flrst order ordinary differen- 
tial equations which can easily be integrated via elliptical 
integrals or Runge-Kutta schemes. The code is very fast 
because the position on the disk (r, 6, <&) is folded di- 
rectly onto the observers screen (x, y) without following 
the path stepwise between these two 2D objects (see Fig. 
0)). This is an alternative to the first historical method 
using t ransfer functions (.Cunningham 1975il wh ich is still 
m use l)Bromlev et al.lll997t (Speith et al.lll995(l or direct 
shooting, a method that follows the complete ray from 
locus of emission to observer. The constraint of dealing 
only with flat disks that lay infinitely thin in the equa- 
torial plane is a good approximation to the vertical col- 
lapsed geometrically thin standard disks, however it is a 
restriction of our code. One recent issue shows the direc- 
tion of future work, where the ray tracing code has to 
be coupled to hydrodynamical or MHD accretion mod- 



6 



A. Miiller & M. Camenzind: Relativistic emission lines from accreting black holes 



els ijArmitage fc Revnoldsll20o'3|) . This coupling allows to 
study the variable emission near black holes. 
The KBHRT code encapsulates all relevant tools to vi- 
sualize relativistic disks and emission lines. One can study 
the distribution of the generalized Doppler factor allover 
the disk or the emission itself. 
The generalized Doppler factor is defined by 

1 



9 '■■ 



I'obs 



1 + 2 



Pobs 
Pcm 



(10) 



where fobs and correspond to observed and emitted 
photon momenta, z to the redshift and p denotes the 
rest frame of the plasma. The Zero Angular Momentum 
Observer (ZAMO) or Bardeen observer momenta have to 
be boosted into the rest frame of the plasma to derive g 
by means of a local Lorentz transformation, 



(t) 



(7 — 1)pq')W^-'^w^'^ — 7p 



(11) 



with the Lorentz factor 7 as measured by ZAMOs, the 
Carter momenta p'^^ and the modulus of the velocity v. 
Only the component p* is needed and one gets 
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.(12) 



With — > 1 for r — !■ cx) as valid for a distant observer, 
the generalized Doppler factor in the emitter frame is eval- 
uated to 

9 = 



7 



(1 - LuX) - 



p 



7 



1 



pVA P 



(13) 



Here, the already introduced Boyer-Lindquist functions 
are used again. Furthermore, there is the plasma veloc- 
ity field in the ZAMO-frame (denoted by round brackets) 
^(r)^^(f) gj^^ two polynomials of the fourth order, 7?,o 
and 9, associated with the Carter Constant and the spe- 
cific angular momentum A = J/E^ the ratio of the two 
conserved quantities along each photon path, angular mo- 
mentum J and total energy E. The polynomial TZ (in most 
general form, see Eq. (|20|l ') becomes for photons (m — 0) 



= r"* + (a^ - - C)r^ + 2 [C + (A - af] r - a^C. (14) 



Usually, one assumes Keplerian rotation for r > rn 
sets 



and 



(15) 



The angular frequency, f2, has for radii greater than 
marginal stability a Keplerian profile 



n = nK = ± 



'M 



for r > r,, 



(16) 





Fig. 5. Simulated disk image at inclination of 30° and 
extreme Kerr a = 0.999999 with color-coded distribution 
of generalized Doppler factor g (white corresponds to 5 = 
1). The beamed blue approaching left side is clearly seen 
(triangular segment, 17 > 1). Gravitational redshift (dark 
colors, g < V) darkens the inner edge until g vanishes at 
the horizon. 

(with upper sign for prograde, lower sign for retrograde 
orbits) and usually one forgets about the contributions 
of w^'-' and v''^\ These components ca n nev ertheless be 
important for the emission line profile l|Miiller 2000). A 
more realistic and complicated velocity field than only 
Keplerian rotation will be discussed in the next section. 
For r < r„is, constant specific angular momentum is as- 
sumed. This is motiyated b y ID hydrodynamic accretion 
disk models llCra,ciall2nni a nd rela tivistic 2D hydrody- 
namic accretion disk models ()Spindeldreher..2002.) . Here, 
one can set 



— rJin — uj 



a Xnis c 

r — tor r < rms, 

1 - wAms 



(17) 



with specific angular momentum at the marginally stable 
orbit given by 



Ams — 



'^ms (^K,ms " ^ms) 



,ms 



(18) 



(The index ms refers to the quantity to be taken at the 
radius of marginal stability, rmg.). 

In Fig. |31 the simulated distribution of the Doppler fac- 
tor g from Eq. (|13|l on a disk with parameters typical 
for Seyfert-1 galaxies (i = 30°) is shown. First, one dis- 
covers the usual separation in blueshifted and redshifted 
segments due to rotation of the disk. The left-to-right 
symmetry is compared to the Newtonian case distorted. 
The left side of the disk rotates towards the observer and 
this radiation is beamed due to special relativistic effects: 
the matter rotation velocity becomes comparable to the 
speed of light. This fact can be derived from Fig. |H1 illus- 
trating the prograde and retrograde velocity component. 
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Fig. 6. Radial dependence of the velocity component u^*-' 
in the equatorial plane for Kerr parameter a — 0.8. 
Prograde case is shown in the upper curve; retrograde case 
is illustrated in the lower curve. The matter starts in free- 
fall from a drift radius of Rt — 3.0 Tg. u*^*) vanishes self- 
consistently (for ZAMOs) at the horizon at rn = 1.6 Tg. 



in the ZAMO frame. The Keplerian profile of w*^*) 
holds only until reaching the radius of marginal stability, 
r-ms- For smaller radii, one has to model this velocity com- 
ponent. The simplest way is to assume free-falling matter. 
In Fig. the drift radius starts already at Rt = 3.0 Tg. 
Then the test particle falls on geodesies onto the black 
hole. The velocity field is a complicated superposition of 
t;^*) and w*^''^. Besides, the plot depicts that the compo- 
nent w*^*) as observed by the ZAMO reaches a magnitude 
that is comparable to the speed of light: symmetrically in 
both cases, prograde and retrograde rotation, the maxi- 
mum rotation velocity is about 0.4c! Additionally, in both 
cases, at the horizon, this velocity component drops down 
to zero. 

For completeness, the distribution of g down to the event 
horizon of the black hole is calculated to illustrate mainly 
the effect of gravitational redshift. Certainly, a disk that 
touches the horizon is unphysical and would not be stable, 
unless in the case of maximum Kerr, a — 1.0, where all 
characteristic black hole radii coincide to 1 rg. Near the 
horizon the g-factor (Eq. (|13|l ') significantly drops down 
and vanishes at the horizon itself (or equivalently, the red- 
shift z becomes infinity): the curvature of the black hole 
grasps everything that is near, even light. Fig.[3now shows 
the Doppler factor to the fourth power. In our approach, 
this is an essential ingredient for any emission originat- 
ing from the vicinity of a black hole (ray tracing using 
transfer functions has a weight to the third power after 
all). This weight to the fourth power shows the relevance 
of the generalized Doppler factor. The g-factor signifi- 
cantly decreases inwards and becomes zero at the event 
horizon. The g-factor to the fourth power shows therefore 
a dramatic tininess around the black hole. T his feature is 
called the central shadow of the black hole ((Falcke et alJ 
boOO^') . Likewise, one remarks a very bright region on the 




Fig. 7. Distribution of the Doppler factor g"^ over the disk. 
The inner disk edge touches at rjn = 1.0015 rg the hori- 
zon of the black hole with a — 0.999999. The inclination 
angle is 30° and the outer edge is at rout — 30.0 rg. The 
shadow of the black hole is clearly seen. Beaming on the 
left approaching side brightens any emission originating 
there. The right receeding part of the disk is darkened. 
This distribution is always folded in the flux integral and 
therefore important for any emission near black holes. 

approaching part of the disk (beaming). For midly and 
highly inclined disks there is a chance to measure this clear 
brightness step, e.g. with sub-mm VLBI in the Galactic 
Center. This may be a method to prove the e xistence of the 
supermassive black hole with 3.7 x IO^Mq llSchodel et alJ 
200^ associated with the source Sgr A* by its radio emis- 
sion. In near future, direct imaging of black holes may 
be possible! 

4. The velocity field of the emitters 

The generalized Doppler factor g depends on the plasma 
motions (here written in the ZAMO frame), g = 
(7(u^''\ w'^®^ w'^*^). The usual issue forgets about the ra- 
dial and poloidal components and derives line profiles 
from emitting plasma with only Keplerian motion, i.e. 
pure rotation. Recent n on-radiative 3D-MHD p seudo- 
Newtonian simulations i Balbus fc Hawlevl l2002l) again 
confirm this inflow and, additionally, the poloidal motion 
of disk winds. Therefore a more complicated plasma kine- 
matics near the event horizon of Kerr black holes is in- 
vestigated in this section. The behaviour of the velocity 
components v^''\ v^^^ and v^^^ at the horizon is presented 
in detail below. 

The radial velocity component will be determined by the 
accretion process. In this paper, it is considered that 
the radial drift starts at a certain radius, the trunca- 
tion radius Rt, that can even be greater than the ra- 
dius for marginal stability r^g. Truncation is an essen- 
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radius [rg] 

Fig. 8. Distribution of specific angular momentum, A(r), 
for Kerr parameter a — 0.8. Fixing a constant specific 
angular momentum, At , is only allowed inbetween the two 
levels: the lower one, Ams, and the upper one, Amb- 

tial feature of radiative accretion theory. This is moti- 
vated by radiative hydrodyn amical accretion models of 
ijHuieirat fc Cam enzind"200rfl . discussed in Sect. |21 The 
correct free-fall behaviour in the Kerr geometry starting 
at radius Rt can be studied by introducing the radial ve- 
locity component 



dr dr dr 
dt dr dt 



n 



A 



.;A 



TP- 1 



;A' 



(19) 



where TZ denotes the radial potential l)Carteilll968jl . This 
polynomial generally depends on Carter's constant C, spe- 
cific angular momentum A, total energy E and mass of the 
plasma test particle m, 



-aA]^-A 



C ,o ^m? 



.(20) 



Carter's constant C vanishes in the equatorial plane, 
u'' is in a normalized ZAMO frame (indicated by round 
brackets) 



exp(/^r) 



a" 



(21) 



This quantity has to be evaluated in the equatorial plane, 
9 = ^. The specific angular momentum, X — J / E, follows 
from 



E 
m 

J 



Vr^ - 3Mr ± 2a\/Mr' 



-•\/r 



(22) 
(23) 



(upper sign prograde, lower sign retrograde orbits). 
This provides the radial drift 



E(l - wA)' 
The kinematics at the event horizon. 



(24) 



(25) 
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Fig. 9. Radial distribution of prograde (solid curve), ^1+, 
and retrograde (dashed curve) angular velocity limit, f2_, 
for Kerr parameter a = 0.8. The dotted curve illustrates 
the frame-dragging potential, lu. At the horizon, here at 
th = 1.6 Tg, all curves coincide. This is the visualization 
of the frame-dragging effect. 

can be explored by evaluating the Boyer-Lindquist func- 
tions 



A(rH) = 0, a(rH) = 0, 

^7,/2{rii) = 2rH, 0J7T/2{rii) = 2, 
^(?'h) = a/2rH = fla- 
This finally yields 

1 



(26) 



2ME~^J) 
1 



E{l + 2^/l-a^ + l)-aJ 



2niE - aJ I 
1. 



2£;(1 + ^/T~^) - aJ 



(27) 



Therefore, matter in the equatorial plane enters the hori- 
zon with the speed of light. This is a general feature of 
the Kerr geometry and fully independent of any accretion 
model. 

A realistic approach to model the radial drift is to fix the 
conservatives E and J when reaching the truncation ra- 
dius Rt with the values that they take at this radius. Then, 
the particle will fall freely into the black hole. The disk 
is truncated at Rt > r^s, the truncation radius gives the 
conservatives, E{Rt) and J(i?t)- Otherwise, the radius of 
marginal stability fixes the conservatives to i?(rms) and 

Jirms)- 

It is not possible to give arbitrary values for constant At , 
that means to fix arbitrary Rt. As ha s been discussed in 
torus solutions on the Kerr geometry ijAbramowicz et alJ 
Il978|) . there is only an allowed interval of Ams < At < Amb 
for the specific angular momentum, where ms denotes the 
orbit of marginal stability and mb the marginally bound 
orbit. This allowed interval At € [A,ns, Amb] is depicted in 
Fig. |Slfor a Kerr parameter of a = 0.8. The general dis- 
tribution of the specific angular momentum in the Kerr 
geometry is represented by the curve. The drift radius has 
to be chosen so that its according specific angular mo- 
mentum. At, just lays inbetween the two horizontal lines. 
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Fig. 10. Radial velocity component in the ZAMO for dif- 
ferent radial drift models, that means different rotation 
state of black hole, a, and drift radius, Bt- In all cases, 
the radial drift becomes light speed at the horizon radii. 
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Fig. 11. Illustration of different radial emissivity models 



These lines belong to the specific angular momenta at the 
radius of marginal stability, Ams, and at the marginally 
bound radius, Amb- The physical reason for this restric- 
tion is, that satisfies the condition 



<n<n. 



a 



il± — UJ ± — 



(28) 



These two limits are shown in Fig. for a = 0.8. It must 
be considered in modelling the radial drift with constant 
specific angular momentum that the resulting fl lays in- 
between these two curves, fi+ and Q^. Otherwise, the 
model is invalid. This is because a global observer is glob- 
ally time-like, UaU" < 0. By the frame-dragging effect, 
this observer becomes light-like, because the horizon is a 
light-like surface {null surface). If one chooses arbitrary 
specific angular momenta, one can run into trouble and 
injure the above allowed stripe for fl. In general, a slow 
rotating Kerr black hole, e.g. low Kerr parameter a, allows 
higher truncation radii. This can be investigated in Fig. 
llUl where models with different radial drift are compared. 
The models differ in Kerr parameter, a, and drift radius, 
i?t, where the radial velocity component, u'^''' , starts to be- 
come finite. As can be seen, the radial velocity increases 
very steeply to smaller radii and becomes always the speed 
of light as has been proven analytically in Eq. H27|l . This 
component is the additive that changes the gi-factor in 
Eq. significantly. The essential consequence is that 
the gravitational redshift effect is enforced what can be 
directly investigated in Fig. El the inner edge shows a 
clear darkening due to smaller g-factors. 
Let us now consider the toroidal velocity component, f . 
From the Carter momenta it follows that 



V 



= -{Q-uj). 



(29) 



Another general feature of the Kerr metric is the frame- 
dragging effect. Anything must co-rotate with the black 
holes horizon when reaching rjj 



(30) 



This can be immediately deduced from the general for- 
mula valid in the Kerr metric 



A 



uj^ 1 — uiX ' 



(31) 



and the fact that a(rH) = 0. So, the azimuthal velocity 
in the ZAMO frame vanishes in general at the horizon as 
can be seen from Eq. H29|l . Fig. confirms this for a — 0.8 
and drift radius _Rt = 3.0 

The poloidal velocity component, v®, in the ZAMO frame 



d9 _ dO dr 
dt dr dt 



p2 1-ujX 



(32) 



vanishes in the equatorial plane, hence the polynomial O 
defined by 



C_ 



"^ ,1 ,? 9 
- 1 I + A^cosec^ 



cos^6' 



(33) 



vanishes for 9 = 7r/2. Therefore, the poloidal velocity com- 
ponent can be neglected for cold, infinitely thin standard 
disks. It may have relevance for slim or thick disks, when 
calculating their line flux. However this requires volume 
ray tracing, covariant radiation transfer and more physics 
as for standard disks what lies beyond the scope of this 
paper. 

5. Emissivity models 

In general, the disk emissivity is the quantity that has to 
be folded into the flux integral to evaluate the observed 
flux from accretion disks. One applies the emissivity as 
kind of profile function that has a radial and potentially 
an angular dependence. The resulting spectra (continuum 
flux, line flux) have therefore no physical flux units but 
arbitrary units. 

The classical approach to model emissivities of ac- 
cretion disks arose in the r elativistic generaliza- 
tion of (non-relativistic) SSDs ((Page fc Thornd Il974l: 
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iNovikov fc Thorndll974 . Here, the radial emissivity pro- 
file follows a single power law with slope index j3 = 3.0 

e(r) cx r"". (34) 

This was motivated by the physics of standard accre- 
tion disks where other physical quantities follow this 
power law behaviour, too. Classical standard disks trun- 
cate at the radius of marginal stability, Tms, where con- 
sequently the emission of the disk breaks off naturally. 
In the historical issue l|Page fc Thornelll974 the Zero- 
torque boundary condition (ZTBC) has been applied. 
Then, the emissivity decreases steeply at this charactcr- 
istic radius. Modern acc retion theory resigns the ZTBC 
l)Merloni fc Fabianllioosfl . On the one hand, this was the- 
oretically motivated by efficient mechanisms to produce 
MHD turbulen ce, the well-known magne to-rotational in- 
stability CMRI. iBalbus fc HawlevI l|l99l|) ). The MRI and 
the coupling of magnetized accretion disk to the black hole 
magnetosphere induces magnetically driven torques that 
injure the ZTBC. On the other hand, observations sug- 
gest significantly deviating disk emissivity profiles, e(r) cx 
r~^'^, giving a hint to higher dissipation at radii r < r^s 
llWilms et al.ll2n0ll) . 

The classical ansatz, /3 = 3, was generalized later to a 
single power law index different from 3.0. The index (3 
may vary locally on the standard disk because different 
illumination from the hot corona and complex structure 
of the ionization layer on the disk can cause changes in 
the emissivity: /? becomes a fit parameter in X-ray spec- 
troscopy. Another aspect is, that at r > r^^s the disk may 
be truncated. The inner accretion disk topology depends 
mainly on the total mass accretion rate as investigated in 
hydrodynamical accretion theory. For high accretion rate, 
the disk extends nearly to th e black hole horizon and the 
ADAF (lNa,rava,n fc V il li 994ll is small. For low accretion 
rate, the ADAF is larger: a spheroidal evaporated hot re- 
gion is attached to the standard disk and forms a sandwich 
configuration. A variable accretion rate is responsable for 
the different spectral states (very high, high, intermediate, 
low, quiescent) and radia lly oscillating transition radiu s 
between SSD and ADAF l|Esin et al.lll997t ICraciJ2002t) . 
In the TSD scenario that emerges only in radiative accre- 
tion theory, it is possible that the matter then decouples 
from the standard disk in free-falling clouds. This plasma 
packages drift radially inwards and are advectcd from the 
black hole. This is major motivation of our new plasma ve- 
locity field model presented in Sect. 21 This phenomenon 
justifies an emissivity profile that does not decrease as 
sharp at radii comparable to r^as- 

Furthermore, one can test so-called double power laws 
or broken power laws where at a certain radius the 
power law index jumps to another value. 
In a simple approach, it is easy to modify the classi- 
cal emissivity profile in the first step by a cut— power 
law. The most simple model would envisage a direct 
vertical cut in the emissivity profile. This may be over- 
simplified because it is expected that the emissivity is 
enhanced at the inner truncation due to stronger illu- 



mination by the corona. An even more adequate model 
considers the single power law with a modification by 
an exponential factor that suppresses the emission at 
smaller radii. Then, the emissivity profile decreases more 
smoothly at the truncation radius. Physically, this is mo- 
tivated by pseudo-Newtonian, radiative two-component 
hyd rodynamics accretio n disk models in 2.5D axisymme- 
try ijHuieirat fc Camen zind 2000) where Comptonization, 
synchrotron radiation, bremsstrahlung and ion/electron 
conduction are considered. These simulations have shown 
that disk truncation at radii r > seems to hold. It 
is supposed that this scenario is realized in AGN, miqro- 
quasars and GBHCs. 

Then, the radial emissivity in cut-power law form satisfies 



e(r) cx exp 



Rt 



(35) 



with truncation radius Rt , steepness parameter a and sin- 
gle power law index f3. The maximum of this distribution 
can be found at r^ax = RtOi/ (3. If one assumes a standard 
emissivity power law, (3 = 3.0, for the second factor, one 
can fix a = 3.0, too. Then it is ensured that the emissivity 
peaks at Rt and exponentially decreases at smaller radii. 
This disk emissivity profile is comparable to the historical 
work ijPaee fc Thornell974|) but exhibits a higher emissiv- 
ity at the inner disk edge. This is plausible because X-ray 
irradiation may be enhanced due to better illumination 
for geometrical reasons. 

But there are alternatives to model the emissivity pro- 
file. In a new approach, a localized emissivity realized by 
a Gaussian emissivity profile is assumed 



e(r) cx exp 



(r-i?t 



(36) 



One can easily control the width of the Gaussian by the 
parameter a-^. Narrow and broad Gaussians with constant 
(Tr are investigated. Another possibility is variable width 
by assuming a-^ = rjRt. Here the width scales with the 
truncation radius Rt, rj denotes a constant factor that can 
be suitably chosen. 

Gaussian emissivity profiles serve for modelling emission 
from a ring. In principle, the ring is the Greens function of 
the axisymmetric emission problem in the Kerr geometry. 
One can simply imagine that any spectral shapes origi- 
nating from extended axisymmetric and fiat emitters are 
composed of small rings, e.g. infinitely narrow Gaussians. 
Hence, one can produce a relativistic emission line profile 
by superimposing the fiux data of several rings. Therefore, 
the Gaussian emissivity profile is well-suited to study dif- 
ferent emitting regions. 

Fig. ^2 illustrates all radial emissivities presented here in 
direct comparison. 

6. Visualization of standard disks 

With these preparations, the investigation of the emission 
distributed over the disk is straight-forward. One very in- 
teresting case is to study the influence of the radial drift 
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Fig. 12. Influence of radial drift. Top: pure Keplerian ro- 
tation w^*) ^ 0, v^'^^ — 0. Bottom: Keplerian rotation plus 
radial drift u*^*^ ^ 0, v*-"^^ ^ 0. Enhanced gravitional 
redshift is clearly seen by darkening near the horizon. 
Parameters were a = 0.1, i = 40°, Hn = = 1.996 rg, 
''out = 10-0 Tg. The radial drift starts at i?t — 5.0 Vg 

i'^''-' ^ 0, as presented in Sect.^ Eq. I|24|) . This new ansatz 
to model the velocity field results mainly in an enhance- 
ment of the gravitational redshift process (compare Fig. 

Another instructive example are high— inclined disks. 
Surely, this can not be applied to AGN type-1 and AGN 
type-2 are obscured by the large-scale dust torus, so that 
there is few chance to observe broad X-ray emission lines 
from type-2 AGN. But maybe this can be applied to stel- 
lar black holes in suitable orientation that can be found 
in microquasars and GBHCs. Besides, these studies are 
interesting to understand GR effects in this part of the 
parameter space. 

A narrow ring laying directly around a Kerr black hole 
horizon is examined. Fig. 1131 now shows firstly the distri- 
bution of the generalized Doppler factor g from Eq. (|13|l 
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Fig. 13. The upper disk shows the distribution of the g- 
factor over a narrow ring with radial drift starting at Rt — 
3.0 rg. The inclination angle is 80°, the black hole rotation 
is fast, a = 0.8, inner edge is the horizon at Hn — rn — 
1.6 Tg and outer edge is only at Tout = 4.6 rg. The lower 
disk indicates the associated distribution of the emission of 
the same narrow ring. Here, the emissivity follows a broad 
Gaussian profile with ~ 3.0 rg and constant Gaussian 
width cr,. — 3.0. 

under consideration of a Keplerian velocity component, 
i;*-*-*, and a drift component, v^'^^ that is switched on at 
truncation radius Rt = 3.0 rg. Generally, for highly in- 
clined disks the beaming of the approaching part is very 
strong. This illustrates the deep blue spot at the left disk 
segment in Fig. E| The left approaching part of the disk 
exhibits a clear brightness step as can be investigated in 
the emission distributed over the disk. Apparently, this is 
a consequence of the g-factor jumping from high values in 
the beaming feature to very small values at gravitational 
redshift feature respective the horizon. If observers may 
find a possibility to look at high-inclined standard disks 
of cosmic black hole candidates, this step would be a dis- 
tinct observable feature confirming the Kerr geometry^. 

7. Relativistic emission lines: calculation and basic 
features 

Let us start with the calculation of the line fiux from the 
rendered disk images. Traditionally in such calculations, 
it has been assumed that the emitting volumes are embed- 

^ More images of disk distributions and 
emission line profiles can be found at 
http://www.lsw.uni-heidelberg.de/users/amueller, 
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ded into the surface of the disk and determine the gener- 
ahzed Doppler factor for each pixel in the image. So, in 
principle, one has to integrate over the disk images (that 
show the emission on each pixel), because a distant ob- 
server is not able to resolve the emitting disk. Physically, 
the observed spectral flux is evaluated generally by 



Beaming 



F' 



obs 



'— ' robs 



dE i: 



(37) 



with the observed intensity 7°^" and an element of solid 
angle dS, covering a pixel of the disk as seen by a distant 
observer. The observed intensity into the rest frame of the 
emitting plasma (denoted by ~ ) holds 



3 rem 



(38) 



with the generalized Doppler factor, g, as introduced in 
Eq. H13|l . In the emitter frame, one has to make an ansatz 
for the intrinsic line shape. Usually, a Dirac delta func- 
tion suffices to describe the transition at rest frame fre- 
quency i/Q. One can assume a Gaussian profile motivated 
by Compton broadening, too. This is only one part, the 
frequency dependence. The other part is the radial disk 
emissivity that satisfies one of the models elaborated in 
Sect. El In our work it is consecutively assumed that 



eir)S (i^om - fo) 



(39) 



Finally, one has to consider a special feature of the delta 
distribution 

S{vcm - I'd) = 5{{uohs - gi^o)/g) = gSivohs - g'^o)- (40) 

So, plug all this into the flux integral (|37|) one gets: 



Fobs {Eohs 



(41) 



image 



with observed flux, i^obs, line emission energy in the local 
rest frame, i?o, observed line energy in the observers frame, 
i?obs, and solid angle element, c?S. Eq. (|41() shows that the 
flux is not a trivial integral over the disk but a complicated 
convolution of intrinsically emitted line weighted with the 
generalized Doppler factor to the fourth power and the 
radial emissivity function. 

Line profiles are simulated that depend in general on the 
parameter set {a, z, rjn, 7'out, i''-*-' , w'^''^ w^®^ ,e{r)}^ with the 
Kerr parameter a, the inclination angle of the disk i, the 
inner Hn and outer edge rout of the thin standard disk, the 
plasma velocity field in the ZAMO frame , w^"'-' , v'^^^ as 
elaborated in Sect. Eland the emissivity law e(r) that can 
depend on one, two or three parameters itself, as has been 
shown in Sect.El 

Therefore, one deals with a rather huge parameter 
space that results in a variety of emission line shapes. 
Consideration of astro-chemistry makes it even more dif- 
ficult, because different species and transitions (Fe Ka, Fe 
K/3, Ni Ka, Cr Ka etc.) can contribute to the observed 

X-ray line feature at approximately 6.5 keV. 

Generally, there are three effects I Fabian et alJll98^ 




Doppler 
Horns 



Gravitational 
Redshift 



Fig. 14. Three effects form the line profile. This simulated 
profile of an iron K line has typical parameters for Seyfert- 
Is. The observer is intermediately inclined, the plasma 
rotates only Keplerian and typical disk sizes fin = Tmsj 
''out — 50 r„ were taken. 



that influence the emission line profile^: the Doppler 
efTect, already known from classical Newtonian emit- 
ting disks, forms a characteristic double horn structure 
at medium to high inclinations. At low inclinations, the 
profile is triangular. The Doppler feature is for relativistic 
accretion disks distorted. 

Beaming is a special relativistic effect. The rotating mat- 
ter moves relativistically fast and the radiation is coUi- 
mated (beamed) in direction of motion. Beaming intensi- 
fies the blue wing of emission lines with respect to the red 
wing. This effect is extraordinary important for disks with 
inner edges close to the event horizon because there the 
toroidal velocity as seen by ZAMOs is comparable to the 
speed of light. At the horizon itself, the toroidal velocity 
steeply decreases to become zero at th as shown in Fig. 
Gravitational redshift is a general relativistic effect. 
Matter and radiation near the black hole feels the strong 
curvature of space-time and looses energy when escaping 
this highly-curved region. The consequence is a shift of all 
photons that succeed in escaping the black holes sphere of 
action to the red branch of the spectrum. This effect re- 
sults in an elongated red tail of the line. Besides, it lowers 
the relic red Doppler peak with respect to the blue wing. 
The frame-dragging effect in the Kerr space-time forces 
anything to co-rotate. In particular, at the rotating hori- 
zon anything co-rotates with the angular velocity fin — ^ 
dictated by the black hole, compare Eq. (|30|l . This means, 
that photons that do not rotate in the black holes direc- 
tion may be forced to turn back. In Schwarzschild space- 
times, a — Q, rotating matter is forced to stop rotation at 
the horizon, building up a boundary layer. 
The combination of all these effects deform the spectral 
emission line to a very skewed profile with typically a long 
red tail (gravitational redshift), an intense beamed blue 



Thermal broadening has a tiny contribution and can be 
neglected. 
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peak (beaming) and a double-peaked structure (Doppler 
boosting) as illustrated in Fig. But the topology of 
the line depends strongly on the inclination angle and the 
emissivity law (see Sect-EJ. 

The parameters used to fit an relativistic emission line 
profile are numerous. In general, the strongest influence 
on both, broadness and topology, has the inclination of 
the disk. At low inclinations, the observer does not feel 
anything of the rotation of the disk and stares onto a 
nearly face-on oriented disk. There is no significant ve- 
locity component towards the observer. The emission line 
profile is triangular. This changes dramatically at high in- 
clinations: the observer sees the approaching and receed- 
ing part of the rotating accretion disk. One can easily fix 
the inclination in a first approximation by determining the 
high-energy cut-off at the blue wing of the observed line. 
The emission line profile is double-horned. 
Another aspect is the rotation of the black hole. For phys- 
ical reasons, it seems clear that the Schwarzschild solution 
is ruled out. From stellar evolution theory it is known that 
stellar black holes formed by gravitational collapse of a 
rotating progenitor stars. Angular momentum can not ra- 
diated away by emission of gravitational waves. Therefore, 
at least stellar black holes are supposed to rotate. Even 
though the formation process of supermassive black holes 
is still enigmatic, one can assume that they form by merg- 
ing of stellar Kerr black holes (hierarchical growth) or 
other rotating constituents. Then, it is very likely that 
black holes in general carry angular momentum. And even 
if black holes do not rotate from the beginning of their 
formation, the accretion mechanism makes sure that the 
Kerr solution is the approbiate metric: matter carries an- 
gular momentum and winds up the black hole spin when 
accreted (compare Sect 13)) ■ 

But the precise rotation state is unsure: it is not clear 
whether to take a Kerr parameter of 0.5 or 0.998. It may 
be a good approach to assume high rotation near extreme 
Kerr (a w 0.99) because these 'elder black holes' accrete 
already for a long time and are wind-up strongly. Besides, 
the relativistic iron lines suggest that the emission orig- 
inates from the innermost regions of Seyfert galaxies. A 
Kerr parameter of a = 0.5 is associated with an orbit of 
marginal stability of 4.23 rg. This may be to far away 
to cause the measured gravitational redshift effects. Disk 
truncation softens this argument as elaborated in Sect. [5] 
As can be seen in the equation of the generalized Doppler 
factor H13(l . the velocity field of the radiation emitting 
matter is very important for the resulting emission line 
profiles. Typically, Keple rian motion (see E g. has 
been taken into account ijFabian et al.ll200o(l . Here, only 
the component u*-*^ is non-vanishing whereas f and 
v^'^^ are zero. But it is certainly essential to incorporate 
a radial drift, e.g. a free-fall motion that is comparable 
to the modulus of the Keplerian motion when reaching 
the orbit of marginal stability. Here the correct fully rel- 
ativistic free-fall formula in the Kerr geometry, Eq. H24() . 
is adapted. 

Poloidal components may also contribute. There is strong 



evidence that the jets of AGN and the "blobs" of GBHCs 
are magnetically driven and form in the direct neighbor- 
hood of the black holes horizon. The interaction of the 
(ergospheric) accretion disk with the black hole magneto- 
sphere drives torsional Alfven waves and finally an out- 
flow (wind) that is coUimated on larger scales to form 
jets. Certainly, the jets of Seyferts are weak and this argu- 
ment holds rather for another type of AGN, the Quasars. 
Besides, to adequately consider poloidal motion, the ray 
tracer has to be modified to a volume ray tracer that ren- 
ders the emission of 3D objects. But this has to be coupled 
to radiative MHD on the Kerr geometry, which is 
still an unresolved problem. 

Recently, the variable iron line profiles from pseudo- 
Newtonian MHD accretion t heory have been calculated 
((Armitage fc R,evnoldsl l2003l) . This model mimics the 
Schwarzschild geometry and shows nicely the depen- 
dence of the line flux from accretion but does still 
not incorporate relativistic effects of a Kerr black hole. 
This issue hints into the right direction of emission 
line simulations: coupling of accretion theory and X-ray 
spectroscopy. Another progress has been done to simu- 
late non-radiative accretion flows on the Kerr geometry 
l|De Villiers fc Hawlevl2003|) . The next step will be to con- 
nect Relativity and radiation transfer to study radiatively 
cooled accretion flows on the background of the Kerr ge- 
ometry. Here, in our first investigation, poloidal velocity 
components of the plasma are neglected. Firstly, the influ- 
ence of radial drift and alternative emissivity models are 
studied. 

One further important aspect is the disk size, especially 
the inner edge. The outer edge is rather irrelevant fo r clas- 
sical emissivities, /? > 2.0 (" Revnolds fc Nowakll2003|) . The 
inner and outer radius of the disk are fixed and afterwards 
the image of this disk is determined. The historical ap- 
proach is based on the famous SSD model and considers a 
thin and cold standard disk that extends outwards to sev- 
eral hundreds gravitational radii (~ 10^ rg) and inwards 
to the marginally stable orbit, r^s- For a rapidly rotat- 
ing black hole (a = 0.998), r^s — 1-23 rg is very close 
to the horizon ru — 1.063 rg! Then gravitational redshift 
influences strongly the radiation from the inner disk edge. 
The accretion disk is supposed to be not as close to the 
black hole. This is motivated by pseudo-Newtonian ra- 
diative hydrodynamics simulations that show for typical 
accretion rates of Seyfert galaxies, M — 0.01 Mq/jx, that 
the disk is truncated at r > r^^- Typical truncation radii 
i?t are 10 to 15 rg, depending on the Kerr parameter a and 
the allowed stripe introduced in Fig. Inland El 



8. Line zoo: topological classification and 
parameter space 

In this section, all types of line topologies emerging in our 
simulations are presented. From the plain line form, one 
can derive the following nomenclature of relativistic 
emission lines 
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Fig. 15. A selection of topology types of relativistic emis- 
sion lines. From top to bottom: triangular, double-horned, 
double-peaked, bumpy and shoulder-like. Line flux in nor- 
malized arbitrary units is plotted over (7-factor. 



— triangular, 

— double-horned, 

— double-peaked, 

— bumpy, 

— shoulder-like. 

Fig. illustrates all prototypes of this terminology. As 
can be seen from the parameters, the triangular form 
follows from low inclination angles at rather high emissiv- 
ity. This is trivial, because the Doppler effect is suppressed 
for nearly face-on observed disks: the red peak of the line 
vanishes. 

Double— horned line shapes are somewhat standard pro- 
files, because many astrophysical objects exhibit these 
typical form. Everything needed for that is a Keplerian ve- 
locity field, intermediate inclination, i « 30°, and a stan- 
dard single power law to reproduce a line profile with two 
Doppler boosted horns, where the blue one is beamed as 
usual. 

The double— peaked profiles are well-known from the 
Newtonian case of higher-inclined radiating disks. Here, 
the width of the relic Doppler peaks is lower as com- 
pared to double-horned profiles. Mainly, this character- 
istic shape is a consequence of the space-time that is suf- 
ficiently flat, so that the typical red tail from gravitational 
redshift lacks. This can be theoretically reproduced by 
shifting the inner edge of the disk outwards. The simu- 
lations with stepwise shifting show nicely the " motion" of 
the red tail until it decreases as sharp as the blue edge. A 
relatively flat space-time is already reached around 25.0 Tg 
corresponding to 0.25 AU for a typical Seyfert galaxy with 
10® M0. Weakly accreting Seyfert galaxies with distant in- 
ner disk edges may fit this scenario. 

Bumpy profiles are mainly caused by steep emissivity 
profiles. These emissivities cut away the emission from 
outer disk regions, especially the beaming segment on the 
disk as depicted in Fig. Therefore, the line profile lacks 
the characteristic sharp blue beaming peak. The recent 
obser vation of Sey fert galaxy MCG-6-30-15 in the low 
state l|Wilms et al.lt2001.') serves as an example of bumpy 
shapes. 

Shoulder— like profiles exhibit a typically curved red 
wing. This feature is very sensitive to the parameters cho- 
sen. Our simulations showed, that only relatively narrow 
Gaussian, broken or cut-cmissivities could produce such a 
distinct feature. The motivation to produce red shoulders 
arose from ano ther observation o f the Seyfert-1 galaxy 
MCG-6-30-15 iFabia n et al.'200^, now in the high state, 
showing such an extraordinary line shape. Blue shoul- 
ders can be produced by significant outflows, e.e 
vanishing poloidal plasma velocity components 
E'oOO'l . Small shoulder- like features attached to the blue 
wing are typically interpre ted as other specie s different 
from Fe Kg s uch as Fe K/3 jFabian et al.ll2002 ^ or Ni Ka 
IjWang et al These transitions have rest frame en- 

ergies above 7 keV. 

Frequently, X-ray astronomers observe narrow but appar- 
ently separated lines that are superimposed to broad line 
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profiles. The interpretation often given is that the two 
components originate from largely separated region: the 
broad line forms in the vicinity of the black hole whereas 
the narrow one is a reflection of X-rays (coming from 
the central engine) at the large dust torus located on the 
kpc scale in AGN. Indeed, the dust harbors the relevant 
species Fe, Ni, etc. Therefore, these features should not 
be observed at microquasars (and stellar black holes in 
general) in default of a dust torus configuration. In simu- 
lations, it is possible to produce these narrow peak com- 
posites by adding an additional narrow Gaussian line pro- 
file without relativistic broadening. 

After this topological classification, one can start to ana- 
lyze the line characteristics numerically. Fig . II 61 illustrates 
a few line criteria that can be proposed to fix the line 
shape features. The utilisation of these criteria may sim- 
plify the comparison of observed and theoretically derived 
emission line profiles. Taking the absissa as a guideline, 
the minimum energy of the line, -Emin, the energy of the 
red relic Doppler peak (if visible), E^p, the energy of the 
blue relic Doppler peak (if visible), -Bbp, and finally the 
maximum energy of the line, i?max, are introduced. The 
two Doppler peaks may serve to determine their energetic 
distance, the Doppler Peak Spacing: 

DPS-Sbp-Srp, (42) 
or more generally in units of the generalized Doppler factor 



line flujc F=. 



DPS 



.9bp - ffip- 



(43) 



On the spectral flux axis, F^p fixes the relic red and i^bp 
the relic blue Doppler peak. To infer an unit-independent 
quantity, one can calculate the flux ratio of these two and 
get the Doppler Peak Ratio: 



DPR - Fbp/Frp. 



Additionally, the line reaches a maximum flux, F„ 
may often match i^bp due to the beaming effect. 
Another essential criterion that is even accessible by ob- 
servation is the line width, in principle the total area of 
the emission line weighted by the energy where the line 
peaks, Eq. In many cases, Eq may coincide with E^p be- 
cause beaming determines the maximum flux at the blue 
peak. The line width is evaluated numerically 



£/0 



(45) 



where AE is equal to the numerical spectral resolution of 
the emission line. More generally, one can express this in 
units of the g-factor to be independent of the line energy. 



(46) 



In both 



.90 

with the corresponding g-factor go = E^^ / E\ 
cases, constant spectral resolution AE-^ = ASj+i = ... = 
AE respectively Ag{ = Ag^+i = ... = Ag is assumed. 
Please consider that this line width has to be distinguished 




line Qf^orgv E 

Fig. 16. Illustration of line criteria to determine charac- 
teristics of a typical emission line profile. The position of 
the relic Doppler peaks (if existing) can be used to fix their 
flux ratio (Doppler Peak Ratio, DPR) and energetic dis- 
tance (Doppler Peak Spacing, DPS). Other specific quan- 
tities are the maximal flux, -Fmax, and the minimum and 
maximum energy of a line , -Emin and i?max- These cri- 
teria can be used to describe a line shape, observed or 
theoretically derived, with a few numerical values. 



(44) ^ 



that 




Fig. 17. Relativistic emission line forming directly in front 
of the horizon. The parameters were chosen to i = 30°, 
Kerr parameter a — 0.998, disk size rin = th = 1.0633 rg, 
''out — 30.0 Tg, radial drift starting at Rt — 1.5 rg and a 
localized emissivity by means of a Gaussian profile with 
Rt = 1.5 Tg and (Ti- = 0.4. This more or less unphysical 
example for line emission should demonstrate the strongly 
suppressed emission due to the shadow of the black hole. 
Even if emission lines could form in the direct vicinity 
of black holes, the high power of the generalized Doppler 
factor g reduces the flux and shifts the line significantly 
to lower energies due to gravitational rcdshift effects. 
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Fig. 18. Investigation of the dependence on the Kerr pa- 
rameter a. The plasma motion is only Keplerian, the in- 
clination is 40°. The black hole spin decreases from top 
to bottom line profile, a = 0.999999, 0.5, 0.1. The inner 
disk edge is equal to the corresponding black hole hori- 
zon, Tin = »'H(a), and the disk edges are chosen so that 
the total disk surface remains constant. 

from the observationally accessible equivalent width. 
Now, only a few selected parameter studies are presented 
because as mentioned before the parameter space is huge. 
The first phenomenon to explore is the frame-dragging 
effect or equivalently the depence of the line shape on 
the Kerr parameter a. Fig. ^] shows only three lines 
at constant inclination angle, i = 40°, and classical emis- 
sivity power law, (3 = 3.0, with variable Kerr parameter 
a = 0.999999 (maximum Kerr), a — 0.5 (intermediate 
case) and a = 0.1 (close to Schwarzschild). The plasma 
kinematics is exclusively determined by Keplerian rota- 
tion. The inner disk edge is always coupled to the ac- 
cording horizon radius (Eq. I|25(l ) and the outer edge is 
chosen so that the net area of the emitting disk is al- 
ways constant. Then, the line profiles are compared and 
one can state that the flux grows with increasing a. This 
is because rin = ru decreases with increasing a and the 
emissivity e(r) cx r^-^ '^ hence increases slightly extending 
to lower radii. This feature is therefore somewhat artifi- 
cial, but may be justified by enhanced dissipation at the 
inner disk edge. In this case, the emissivity profile is sup- 
posed to contribute even for radii r < r^g- One can also 
state a changing line topology in this study: the ratios of 
the relic Doppler peak fluxes, DPR, as introduced in Eq. 
m, DPR(a = 0.1) « 2.0 and DPR(a = 0.999999) « 1.7, 
suggest that beaming is enhanced when the the black hole 
rotates faster. This is plausible by considering the increas- 
ing toroidal velocity, w'*^ respectively 17, with decreasing 
radius. 

In Fig. 1171 demonstrates line dampening by the shadow 
of the black hole. Here the radial drift model is applied, 
too, but a pure Keplerian velocity field does not change 
the following argument. The Gaussian emissivity shape 
provides a suitable tool to easily fix a small emitting ring 
laying around the horizon. The applied Gaussian emis- 




Fig. 19. Reproduction of the line set with variable incli- 
nation, i = 10°, 30°, 60° in Schwarzschild i Fabian et all 
.200(1 . The main difference is that here an additional radial 
drift is adapted whereas lFabian et al.l l|200(1j) only assumed 
Keplerian motion. In both cases the emissivity follows a 
standard single power law. For better comparison, the ab- 
scissa is scaled to the iron Ka line with 6.4 keV rest frame 
energy. 

sivity peaks at i?t = 1-5 and vanishes in principle for 
radii r > 2.0 r^. Therefore, the line form is nearly sym- 
metric. This examples depicts that even if line formation 
could happen within some < Tms, the line feature would 
be strongly suppressed by gravitational redshift. As ana- 
lytically proven in Sect. |3| the (/-factor vanishes at the 
horizon. Therefore, the powers of g (third or fourth, de- 
pending on ray tracing technique) folded into the flux 
integral reduce the observed flux. The consequence is a 
high-redshifted and damped emission line. As can be de- 
duced from Fig. El a Fe Ka line would peak at an ob- 
served energy of approximately only 3.4 keV. The line is 
significantly shifted to lower energies due to gravitational 
redshift. Obviously, it would be hard to observe this weak 
line. But in any case, the fluorescence process is supposed 
to be not possible in this proximity to the black hole be- 
cause the accretion flow is too hot. The line dampening 
due to gravitational redshift is a strong competitive effect 
against mo dels with steep emiss ivity single power laws as 
proposed in lWilms et alJ l)200l|) . 

To demonstrate the influence of radial drift, a line set 
originating from a disk at three different inclination an- 
gles, i — 10°, 30°, 60°, around a Schwarzschild black hole 
that extends from rjn — r^g = 6.0 rg to rout = 30.0 with 
only Keplerian plasma motion is reproduced l)Fabian et alJ 
20001) . The same parameters are adapted but the inner 
disk edge is chosen to be at the horizon. Then, a compo- 
nent u^''^ 7^ 0, e.g. a radial drift that starts at i?t = 6.0 rg is 
considered. The emissivity law is chosen to a single power 
law with /3 — 3.0 without inner cut-off to demonstrate the 
influence of emitting radially drifting plasma for r < r^w 
Fig. 1191 depicts the result. A comparison of both line sets 
confirms that the blue edges remain almost constant. But 
radial drift starting already at the radius of marginal sta- 
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Fig. 20. Relativistic emission lines with radial drift model 
and Gaussian emissivity. The drift radius i?t from the 
plasma velocity field couples to the emissivity. From bot- 
tom to top, the drift radius increases, i?t = 4.0, 5.0, 6.0. 
The black hole rotates slowly, a = 0.1, the inclination an- 
gle is 40°, Tin = rn = 1.995 and rout = 30.0 rg. The 
emissivity is modelled with constant CTr = 3.0. 

bility, This = 6.0 Tg in Schwarzschild, suppresses the red 
relic Doppler peak by enhanced gravitational redshift ef- 
fects. The red wings are in all three cases more smeared 
out with radial drift than in the case where only rotation 
is considered. For observation, this means that radial drift 
produces longer and more intense red tails of the emission 
line. 

Now, the relativistic emission line profiles resulting from 
the new truncated accretion disk model presented 
here is investigated. The radial drift superimposed to the 
Keplerian motion as introduced in Sect. 0] and the new 
emissivity model with Gaussian shapes as has been shown 
in Sect. are implemented. The missing link to this ap- 
proach is the coupling of the Gaussian emissivity to the 
radial drift. This model can be justified by the following 
scenario: radiative accretion disks are truncated by ef- 
ficient radiative cooling. This means that matter at the 
inner edge (at Rt > r,„s) of the standard disk evaporates 
and matter clouds start to fall freely into the hole. This 
can be modelled with the radial velocity component w*-"^' 
as introduced in Eq. (|24|) . At the horizon of the black 
hole, the plasma is mostly advected into the Kerr black 
hole and passes the horizon with the speed of light (an- 
alytically derived in Eq. (|27|l and illustrated in Fig. 1101) . 
The emissivity is modelled in a consistent manner: the 
illumination of the inner edge of the standard disk is ex- 
pected to be increased because the disk exposes a larger 
surface to the primary X-ray source, e.g. a hot inner torus, 
an ADAF or generally speaking a corona. Therefore, the 
emissivity is higher than elsewhere on the accretion disk. 
This can be modelled by Gaussian emissivities that de- 
crease very steeply to both wings of the emissivity shape. 
This is plausible, too, because at smaller radii there is 
only free-falling material that is optically thin and be- 
sides it is too hot to permit the fluorescence transition. At 



larger radii, the emissivity is expected to decrease, too, 
because the the surface of the disk becomes too cold for 
fluorescence. Hence, one can motivate to link the Gaussian 
profile of the emissivity to the drift radius. In a first step, 
broad Gaussians with constant width, — 3.0, are ap- 
plied in all three cases depicted in Fig. 1201 The black hole 
rotation state is constant, a = 0.1. The drift radius varies 
from Rt ~ 4.0, 5.0, 6.0 rg and is inserted in both, radial 
drift (Eq. l(2il) and Gaussian emissivity (Eq. JSEJ). Fig. 
1201 now shows that the line flux increases with increas- 
ing drift radius. The calculated line width, in principle 
the area below the line, even grows by a factor of 2, if 
one compares Rt = 4.0 rg with Rt = 6.0 rg. This is be- 
cause the Gaussian width remains constant whereas the 
emitting ring (limited by Gaussian emissivity) grows in 
surface area. But a changing line topology can be derived, 
too. The red shoulder is elongated at Rt = 4.0 rg because 
of enforced gravitational redshift for smaller drift radii 
respectively truncation closer to the black hole. Besides, 
the simulation yields a slight variation in the line form. 
This can be proven with the quantity DPR, the Doppler 
Peak Ratio which is independent from absolute flux val- 
ues. This is the ratio of the flux from the blue Doppler 
peak versus the red Doppler peak. Surely, this can only 
be calculated if the two Doppler peaks exist. These DPRs 
are DPR(i?t = 4.0) = 3.6, DPR(i?t = 5.0) = 4.0 and 
DPR(i?t = 6.0) = 4.3. So, the red shoulder drops down 
with increasing Rt and the blue flux peak increases by the 
growth of the emitting region. 

An alternative to the above issue, is to replace the 
Gaussian emissivity by the cut— emissivity profile as 
introduced in Eq. H35|l . This radial shape ressembles 
to the coronal emissiv ity profiles recently presented 
ilMerloni k FabianI |2003|) . For comparative reasons, all 
other parameters as in the example with Gaussian emis- 
sivity are maintained and only the emissivity is changed 
to a cut-power law with constant a — 3.0 and the same 
truncation radii set, Rt = 4.0, 5.0, 6.0. The behaviour 
at larger radii corresponds to the standard single power 
law cx r"^ *^. The inferred emissivities are shown in Fig. 
1211 The resulting line profiles with these cut-emissivity 
shapes can be investigated in Fig. |221 For the three lines 
itself, the allover topologies are almost identical, even 
DPS remains constant. But the ratio of the relic Doppler 
peaks slightly increases with Rt- DPR(i?t = 4.0) = 4.95, 
DPR(i?t = 5.0) = 5.16 and DPR(i?t = 6.0) = 5.26. This 
effect is mainly caused by the fall of the red wing driven by 
the radial drift. Remarkable is the fact that the blue edge 
stays absolutely constant. Comparing Figs. 1201 and l22l the 
blue edges remain constant. But, it is stressed that at 
the parameters chosen the Gaussian emissivity produces 
shoulder-like line shapes whereas the cut -power law emis- 
sivity produces double-horned profiles. 
The reduction of the red wing flux caused by including ra- 
dial drift is a generic feature. Fig. 1231 confirms this. Here 
two lines are compared where in one case a pure Keplerian 
velocity field applies and in the second case the radial drift 
model is used. This study reveals immediately the influ- 
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Fig. 21. Cut-emissivity with a = (3 = 3.0 for different 
truncation radii. From bottom to top, the drift radius in- 
creases, i?t = 4.0, 5.0, 6.0. With this choice of a and 
/3, the emissivity peaks naturaUy at i?f For r < Rt the 
emissivity decreases exponentially. 
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Fig. 22. Relativistic emission lines with radial drift model 
and cut-power law emissivity. The drift radius i?t from the 
plasma velocity field couples to the emissivity. From bot- 
tom to top, the drift radius increases, i?t = 4.0, 5.0, 6.0. 
The black hole rotates slowly, a = 0.1, the inclination an- 
gle is 40°, Tin = TH = 1.995 rg and Tout = 30.0 Tg. The 
emissivity is modelled with constant a = [3 — 3.0. 



ence of including plasma drift. In this example, the drift 
starts at Rt = 4.5. One can see that the red part of the hne 
drops down when drift is operating. The interpretation is 
that the gravitational redshift is enhanced by the drift, 
compare Eq. (|13|l . and results in a lower red wing flux. 
This statement holds in principle for all simulated lines 
where the drift model is compared to its pure Keplerian 
counterpart. 

Finally, a more complicated line profile with a few species 
in application to a recent observati on of Seyfert-1 gala xy 
MCG-6-30-15 in the high state l|Fabian et al.ll2002|) is 
demonstrated. The characteristic shoulder-like shape of 
this observation can be reproduced theoretically with the 
following plausible parameters: 

— rapidly spinning black hole, a = 0.8, 

— typical Seyfert-1 inclination, i — 40°, 

— inner edge at horizon, nn = rn = 1.6 rg, 

— outer edge at rout = 20.0 rg, 




Fig. 23. Direct confrontation of relativistic emission lines 
with and without radial drift model. A single-power law 
emissivity with /3 = 3.0 is chosen. The black hole spin pa- 
rameter is a = 0.8, the inclination angle is 40°, Hn = rn = 
1.6 rg and rout = 20.0 rg. The pure Keplerian case, top 
line, is modelled by constant specific angular momentum 
for r < rijis = 2.91 rg. The non-Keplerian case, bottom 
line, includes drift starting at i?t — 4.5 rg. 

— standard Keplerian rotation plus 

— radial drift starting at Rt = 4.5 rg, 

— Gaussian emissivity: Rt = 4.5 rg and Cr = 0.4 Rt, 

— multi-species line system consisting in Fe Ka, Fe K/3, 
Ni Ka and Cr Ka 

The relative line strength of the different species are mod- 
elled in an approximation found in Reynolds (1996). Fig. 
1241 now illustrates a line system consisting of iron, nickel 
and chromium. The total observed X-ray fluorescence fea- 
ture around 7 keV is equal to the sum of these contribu- 
tions. The theoretically inferred line FWHM is very high, 
around 1 keV, wh ich has the same ord er as the k, 750 eV 
as can be found in iFab ian et al ]f2002'). These parameters 
fit nicely the observed shoulder-like red wing and even the 
small blue bump beyond the Fe Ka beaming branch. The 
fiux ratio between shoulder- level and maximal flux around 
a factor of 5 is close to the observed value of 4. The other 
elements seem to be vital to model this observed profile. 
However, chromium seems to play a minor role due to 
weak relative strength. 

9. Conclusions 

It has been demonstrated that radial drift, a non- 
vanishing radial velocity component, is a vital ingredient 
in modelling relativistic emission line profiles. Accretion 
theory demands on such more complicated non-Keplerian 
velocity fields. However, our approach is just a phe- 
nomenological model motivated by simulations of 
accreting black hole systems. Some analytical features, 
as the asymptotical behaviour at the black hole horizon 
are independent of these models. One main consequence 
is, that the radial ZAMO velocity becomes the speed of 
light at the horizon whereas the toroidal ZAMO velocity 
vanishes there. Another essential feature in accreting and 
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Fig. 24. Suggestion of a multi-species line system fitting 
the XMM-Newton observation of MCG-6-30-15 (Fabian 
et al., 2002). A line complex consisting of Fe Ka at 6.4 
keV, Fc K0 at 7.06 keV, Ni Ka at 7.48 keV and one weak 
contribution of Cr Ka at 5.41 keV (all energies measured 
in the rest frame) is simulated. 



radiating black hole systems is the steep decrease of the 
generalized Doppler factor causing a shadow of the black 
hole. Hence, any emission near black holes is strongly 
suppressed. This is the competing dominant effect that 
suppresses as well steep emissivity profiles within the 
marginally stable orbit. 

The truncated standard disk (TSD) hypothesis induced 
by radiative pseudo-Newtonian hydrodynamics is a new 
aspect in accretion theory that also influences the models 
for X-ray spectroscopy, especially the formation of broad 
relativistic emission lines. Future fully relativistic and 
radiative MHD is expected to confirm this. Then, a cou- 
pling of 3D ray tracing techniques and plasma dynamics 
are inevitable to explain the morphology and variablity of 
relativistic emission lines, particularly the iron Ka line. 
Besides, one hope is that these demanding simulations 
will enlighten the open question of the corona topology. 
The answer is essential to understand X-ray irradiation 
of thin accretion disks. 

Alternative disk emissivity models, such as cut-power 
laws and Gaussian emissivities can serve as valuable 
tools to study the huge parameter space of the line zoo. 
The localized emissivity studies with Gaussian emissivity 
laws prove that emission from regions neighboring the 
black hole horizon is massively suppressed by the low 
generalized Doppler factor. Additionally, the exponen- 
tially decreasing cut-power law may be a suitable profile 
function to model the disk emissivities with higher 
magnetically driven dissipation at the inner disk edge. 
A first classification of the variety of relativistically 
broadened emission lines succeeds by a topological crite- 
rion: the lines can be divided in families with triangular, 
double-horned, double-peaked, bumpy and shoulder-like 
line form. It turns out that this first step constrains the 
parameter space significantly. 

Unfortunately, it is still not feasible to solve the fully 
covariant radiation transport problem. Therefore, one 



deals always with some restrictions in the accretion 
models: either the approach is fully relativistic, but 
neglects radiation or the model incorporates radiation 
processes, but mimics relativity only to a certain degree 
by pseudo-Newtonian potentials. The future of accretion 
theory points towards radiative 3D— GRMHD. Until 
this will be possible, one has to extract few informations 
from 'simpler' models with limitations done in the past 
(such as pscudo -Newtonian models , non-radiative 
issues, less dimensions etc.) and puzzle this input to an 
understanding of black hole accretion physics. 
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